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Abstract 

We study the behavior of solitary-wave solutions of some generalized 
nonlinear Schrodinger equations with an external potential. The equa- 
tions have the feature that in the absence of the external potential, they 
have solutions describing inertial motions of stable solitary waves. 

We construct solutions of the equations with a non-vanishing external 
potential corresponding to initial conditions close to one of these solitary 
wave solutions and show that, over a large interval of time, they describe 
a solitary wave whose center of mass motion is a solution of Newton's 
equations of motion for a point particle in the given external potential, 
up to small corrections corresponding to radiation damping. 

1 Introduction 

In this paper we study the effective dynamics of solitary wave solutions of a 
class of generalized nonlinear Schrodinger equations with an external potential. 
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These equations have the form 

i^ = (-A + y)^-/W, (i.i) 

where if> : JBL d x R i— > C, x denotes a point in space WL d , t G R is time, <9 t = t|, 
A = ^j=i ^ i s the spatial Laplacian, V(x) is the external potential and is 

a real-valued, bounded slowly varying function on R d , and / is a map from 
the co mplex Sobolev space Hi(R d ,C) to H_i(R d ,C) such that /(0) = 0, and 
/(VO = describing a nonlinear "self-interaction". Precise assumptions on 

and / will be given in Section |21 Examples of nonlinearities of interest include 
local nonlinearities such as 

f(iP) = \\if)\ 2s iP, 0<s<| A>0, (1.2) 

and Hartree-type nonlinearities 

= A(jy*H 2 )V>, a>o, (i.3) 

where W is of positive type, continuous, spherically symmetric potential function 
which tends to 0, as \x\ — > oo and W * g(x) := J W(x — y)g(y) d d y, denotes 
(spatial) convolution. Of course, A can be scaled out by rescaling if). 

We assume that the nonlinearity in (jl.lj) is such that the Cauchy prob- 
lem has a unique, global solution, if>(x,t), in the space C(R + ; Hi(R d , C)) fl 
C X (R + ; H_!(R d , C)), given an initial condition ip(x, 0) = if) {x) G Hi(R d , C). Re- 
sults on the Cauchy problem associated with (jl.lj) can be found in ^HJ I2SJ E] 
(see Section |2 for a discussion). 

For V = 0, eq. (jl.lj) is the usual generalized nonlinear Schrodinger (or 
Hartree) equation. For self-focusing nonlinearities (as in example (jl.2j) or (jl.3j) 
with W of positive type and A > large enough), it can have stable solitary 
wave solutions of the form 

Va {x,t) := e^-^^ix - a), (1.4) 

where a := {a,t>,7,/i}, and a = v t + ao, 7 = fit + ^-t + 70, with 70 G [0, 2tt), 
ao, v G R d and /1 G R + , constant, and where rj^ is a positive solution of the 
nonlinear eigenvalue problem 

(-A + ^-j^HO (1.5) 
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(see Section |2I). Solutions of eq. of the form (jl.4|) describe solitary waves 
traveling through space with a constant velocity v, and with an oscillatory phase 
given by fit — \v 2 t. Existence of such solutions, for a large class of nonlinearities 
has been established in gU H EJ El 123 I2E1 d • See Section H for an outline of 
results relevant for this paper. 

In analyzing soliton-like solutions of eq. we encounter two length scales: 
the size oc /i" 1 / 2 of the support of the function which is determined by our 
choice of initial condition ipo, and the length scale oc (sup | VV^I)" 1 over which 
the external potential V varies appreciably. We will assume that the ratio 

sup|W(aQ| 

£v = — , (1-6) 

is small; i.e., that the potential V varies little over the support of a solitary wave 
solution. 

When V does not vanish, the wave T] a (x) in (jl.4j) does not solve eq. (jl.lj) . 
However, we expect that if the initial condition 

Mx) :=^(x,Q) (1.7) 
for (jl.lj) is close to rj ao (x) for some cr , in the sense that 

lle-^-M-i&JIk <e , (1.8) 

then for all times < t < — 7—5, where T is some positive constant, the solution 
ip{x, t) remains close to a solitary wave of the form (jl.4j) for some time- dependent 
parameters /j,,v,a and 7 determined by V and ipo- We will show, more specifi- 
cally, that these parameters can be chosen to be solutions of the following system 
of ordinary differential equations 

a = v, £ = -W(a), A = 0, (1.9) 

and 

j = fi+^-V(a). (1.10) 

(up to error terms of size 0(e v 2 +e 2 )), with initial conditions given by a(0) = a , 
v(0) = v Q , /i(0) = fi , and 7(0) = 70, with a ,fo,/io;7o as i n P-4j) . We observe 
that the first two equations in (jl.9j) are Newton's equations of motion for the 
trajectory (a(t),v(t) = a(t)) of a point particle of mass | moving in the external 
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potential V(a). The center of a solitary wave solution of follows this 

trajectory, up to deviations 0(e v + e ) due to "radiation damping" . 

We state here rigorously the result discussed above for a special class of 
nonlinearities. The general class of nonlinearities is introduced and discussed in 
the next section. We assume the external potential V(x) satisfies the conditions 

V G C 2 and \d%V(x)\ < C a e v |a| , for \a\ < 2. (1.11) 

In e v > is the small parameter introduced in ()1.6|) . In other words, we 

find it convenient to fix the size of the support (the 'width') of the solitary wave 
solution at 0(1) and assume that the external potential V(x) varies slowly. Let 
e:=e v + e . 

Theorem 1. Assume that the nonlinearity f is given by (|1.2|) . and assume that 
the potential V satisfies (|1.11|) with e v < 1. Let I be any closed, bounded 
interval in (0, oo). Let e <C 1 and the initial condition ip satisfy 

lle-^^o-^cJik <eo, (1-12) 

for some uo G M d x ]R d x [0, 2tt) x I . Then there is a constant T > 0, independent 
of e v> e o but possibly dependent on I , such that for times < t < T(e v + e 2 ) -1 , 
the solution to (|l.lj) with this initial condition is of the form 

if,(x,t) = e^ v < x - a)+ ^(ri^x - a) + w(x - a,t)), (1.13) 

where 

\\w\\ Bl = 0(e), (1.14) 
and where the parameters v, a, 7 and fi satisfy the differential equations 

i) = -(W)(a) + 0(e 2 ), (1.15) 

a = v + (D(e 2 ), (1.16) 

7 = A i-V(o) + ii; 2 + 0(c a ), (1.17) 

fi = 0{e 2 ). (1.18) 

The same conclusions hold for nonlinearities of the form f(ip) = gd'ipl 2 )'^ + 
(W * IV'I 2 )^ where W and g satisfy explicit conditions (see the discussion of the 
conditions in Section EJ), provided an additional spectral condition is satisfied 
(see Condition |(F)| in Section |2J). 



4 



The first result of this type was proved by Frohlich, Tsai and Yau "To! IT7] 
for the Hartree equation ( (jl.ip with under a spectral condition (see Con- 

dition^)] of Section 2). The choice of the Hartree type nonlinearity plays an 
important role in [T"""|["[7|. For local, pure power nonlinearities and a small pa- 
rameter, e v , Bronski and Jerrard [7] have shown that if an initial condition is of 
form (jl.4|) . with t = 0, then the solution ip(x,t) of eq. satisfies 

e v - d \ij(^,±)\ 2 d d x ^ IMl^ (1.19) 

in the C 1 * (dual to C 1 ) topology, provided a(t) satisfies the equation = VU (a), 
where V(x) = U(e v x) (see (Q|) ). 

Our approach is built on important developments in the nonlinear Schrodinger 
equation (NLS) in the last 20 years (see [TTJ H3] for reviews). We outline the 
landmark developments briefly here. Orbital stability of NLS solitary waves for 
V = was proved by Cazenave and Lions [12] and M. Weinstein [17J EH] , whose 
result was significantly extended by Grillakis, Shatah and Strauss [20112]. The 
next significant step was made by A. Soffer and M. Weinstein [3B] who proved, 
under some restrictive conditions, asymptotic stability of nonlinear ground states 
for eq. ()l.lj) and by V. Buslaev and G. Perel'man [H] who, motivated by [35] . 
proved asymptotic stability of NLS solitary waves (V=0) in one dimension, again 
under certain restrictive conditions. These results were significantly extended by 
Tsai and Yau [II I?"J iB] , Soffer and Weinstein [S3 ID], Cuccagna [131 EE], Bus- 
laev and Perel'man ||, and Buslaev and C. Sulem "TUj. Furthermore Perel'man 
[35| . and Rodnianski, Schlag and Soffer jHT], have obtained the first results on 
soliton scattering. Many of the issues touched upon in the present paper were 
studied also in work of Gustafson and Sigal [23] on dynamics of magnetic vor- 
tices. 

We also mention interesting non-rigorous results by Pelinovsky, Afanasjev 
and Kivshar [33] and Pelinovsky and Grimshaw [31] on dynamics of NLS solitons 
near the borderline for stability. Nonrigorous results on dynamics of 'center of 
mass' of solitons were obtained by E. van Groesen and F. Mainardi and G. Derks 
and E. van Groesen [2*2*1 E3 • 

Our paper is organized as follows. In Section [21 we present the general hy- 
potheses on the class of nonlinearities and formulate our main result under these 
hypothesis. In Section "31 we explain the Hamiltonian and variational aspects of 
the dynamics given by eq. (jl.l|) which play a role in our proof. In Section [""J 
we describe the symmetries of eq. (11.1}) for V = and find the "zero modes" 
associated with these symmetries. Furthermore, we introduce and analyze a 
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finite-dimensional manifold, M s , of stable solitary wave solutions to eq. (|1.1J1 . In 
Section |5J we introduce a convenient parameterization of functions in a small 
neighborhood of M s in phase space. In Section |BJ we transform the equation 
()l.ip to a moving frame, and then rewrite the resulting equation in terms of the 
parameters introduced in Section El In Sections [7| and |Hl we control solutions of 
our equations of motion in a moving frame by constructing an approximately 
conserved Lyapunov functional. The proof of our main result is completed in 
Section |§1 Some material of technical or review nature is collected in four ap- 
pendices. 

Remark about the notation: in this paper we consider equations, maps and 
functionals on complex spaces such as Hi(R d , C), which sometimes are identified 
with real spaces; e.g., H 1 (M d , R) = Hi(R d ,R) © H 1 (R d , R), under the identifi- 
cation ip <-> (Reip,lmip). In this case the operator of multiplication by i" 1 is 
identified with the operator 



(which defines a complex structure on the corresponding real spaces). Thus a real 
function 77^ is sometimes written as (77^, 0) and similarly an imaginary function 
as 177^ as (0,7/^) (or even as Jr]^). Frechet derivatives are always understood to 
be defined on real spaces. They are denoted by primes, c will denote various 
constants which may depend only on the interval Iq for the parameter /xq. 
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In this section we formulate general assumptions on the nonlinearity / in (jl.lj) 
and state our main result under general assumptions. Our assumptions are 
rather abstract, each responsible for certain aspect of the problem and of our 
approach. Then we discuss specific nonlinearities for which these assumptions 
are satisfied 




(1.20) 



2 Main Result 
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(A) . (Energy) There exists a C 3 -functional F : Hi i— > M such that F'(ip) = f(ijj), 

su P|MIh 1 <mII^'V)IIb(h 1 ,h_ 1 ) < oo and sup| H | Hi < M HF'^^Hh^bch^H-i) < 
oo (here B is the space of bounded linear operators); 

(B) . (Symmetries) F(Tip) = F(ip), where T is a translation 7^ tr : u(x) t— > 

u(x — a) Va G M. d , a rotation 7^ : u(x) i— > u(R~ l x) WR G SO(d), a gauge 
transform : m(x) i— > e 17 u(a;) V7 G [0, 27r), a boost transform Zj 5 : tt(a;) 1— ► 
e5"' 3: «(x) Vf G R d , or a complex conjugation T c : it (a;) 1— > u(x); 

(C) . (Existence of solitons) There is an interval I C K such that V/i G / 

eq. (jl.5j) has a positive, spherically symmetric, L 2 fl C 2 solution 77^, such 
that || |x| 3 ?7 M || + || |x| 2 | V^ M | || + || |a;| 2 (9 M r7 M || < 00 V/i G / (here || ■ || denotes 
the L 2 norm); 

(D) . (Orbital stability) The solutions 77^, /i G /, of described in |(C)[ satisfy 

d,J V ld d x>0; (2.1) 

(E) . (GWP) Eq. (JUU) is globally well-posed in Hi and in H 2 ; 

(F) . (Null space condition) V/x G I, 

N(£ 7? J = span{(0,^), (9^,0), j = (2.2) 

where := —A + /i — /'(t?^), the Frechet derivative of the map t/> h- > 
(-A + ^)V-/(V) at v 

We now discuss conditions |(A)f|(F)| beginning with general remarks. Condi- 
tion |(A)| allows us to define the conserved energy or Hamiltonian 

UvW) = \fw\ 2 + V\^\ 2 ) d d x - F(tp), (2.3) 
and gives us the following estimates on the nonlinearities 

|itf H| < c(M)W| Hl , \R®(w)\ < c{M)W\k> ( 2 - 4 ) 

and 

IIA^IIh^CWHIL (2-5) 
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for any r\ G H 2 (R d ) and w e Hi(R d ). with ||^|| Hl + |M| Hl < M. Here 

R (2) (w) :=F(ri + w)-F(ri) - (F'(r]),w), (2.6) 
R®(w) := F( V + w) - F(r)) - (F'(r)),w) - ^(F"(r])w,w) (2.7) 



v 

and 

Note that 



N v (w) = F'{rj + w)- F'(rj) - F"(rj)w. (2.1 



N v (w)=R^\w). (2.9) 

We can assume without loss of generality that /'(0) = 0. Then one can show 
that J C M + . 

As was mentioned in the introduction, the nonlinearities of interest to us are 
local nonlinearities, 

fmx) = h(\^(x)\ 2 )^(x), (2.10) 
for some real function h on R + , and Hartree-type nonlinearities, 

ffy) = (W*\ij>\ 2 )iJ> 1 (2.11) 

where W is a fixed, real valued, spherically symmetric function, tending to at 
oo. More generally we consider nonlinearities the form 

f(ip)(x) = h(\^(x)\ 2 )^(x) + (W* \ip\ 2 )(x)ij(x). (2.12) 

We discuss now under which conditions on h and W in ()2.12|) Conditions |(A)[|(Fj| 
are satisfied. 

Condition (A) For nonlinearities of type (J2.12)) the functional F is given by 

FW = \j HM 2 ) + \{W * h/f) d d x } (2.13) 

where H(s) = h(p) dp. The functional F is C 3 with the inequalities specified 
in Condition |(A)| satisfied, provided h(s) is C 2 with h^ k \s) < c(l + s q ~ k ) (k = 
0, 1, 2), for some q<2/(d- 2), and d < 4, and W e L p + L°° for some p > d/2. 

Condition |(B)[ This condition is trivially satisfied for ()2.12|) with W(x) = 
W{\x\). 
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Condition |(C)[ Existence of a positive, spherically symmetric solitary wave 
solution to Eq (jl.5j) . was proved in jU E] for local nonlinearities (|2.1(J|) satis- 
fying 

— oo < lim h{r) < fi, 
-oo < lim r~ a h{r) < C, (2.14) 

r— *aa 

for < a < 2/(d - 2), if d > 2 and a G (0, oo), if d = 1, 2, 



3£ > 0, such that / /i(r) dr > 



Additional results, for a large class of nonlinearities can be found in [JTJ H l2~4"l 

121 p. 

For nonlocal nonlinearities (|2.11j) (Hartree equation) with 

W G Lf oc , p > |, W -> as |x| -> oo, (2.15) 

existence of solutions was proved in (TTJ H3 123 I2H U]- 

Moreover, for nonlinearities (|2.10j) or (|2.11|) satisfying conditions (|2.14j) or 
f!2.15|) . solutions to ()1.5|) (solitary waves) are exponentially decaying at oo as 
0( e -VPM); see HE21- 

Condition |(D)| Condition (|2.1|) is a sufficient condition for orbital stability 
of solitary waves for (generalized) nonlinear Schrodinger equations (see [201 IS])- 
This condition is to be checked for each nonlinearity. In the special case of a 
pure power nonlinearity, f{jp) = X^^ip, we have rj^x) = /i2^^ M= i(xy / /I). Thus 
condition |(D)| reduces to the condition that s < 2/d. 

Condition |(E)| Assume that / is of the form (|2~T2|t with h : [0, oo) h-> R, 
smooth and satisfying h(0) = 0, and 

\ti(r)\ < Cil + r**- 1 ) , (2.16) 

for some a G [0, tr^), for d > 3 and a G [0, oo) for d = 1, 2 and 

fr(r) < C(l + r") , (2.17) 



for some (3 G [0, |). Furthermore assume that : R ' i— > R is an even potential 
such that 



G L 9 (R a! ) + L°°(R a! ), (2.18) 
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for some q > 1, q > d/4 and let W + = max(0, W) s.t. 

W+ G L r (R d ) + L°°(R d ), (2.19) 

for some r > 1, r > d/2 for d > 2 and r > 1 for d = 2. Then (jl.lj) is globally 
well-posed in Hi and H 2 (see [TT1 Chapter 6]). 

Condition |(F)| This condition is more delicate. First we observe that 

(9^0),(0 1 ^)eN(A,J ) Vj = l,...d, (2.20) 

due to the fact that r/ M breaks the translation and gauge symmetry of (jl.5|) (see 
Section HI). 

Now we list some facts which are discussed in AppendixO If rj^ is spherically 
symmetric and / is a local nonlinearity, f(ip)(x) = f(ip(x)), then C v can have 
at most one zero eigenfunction in addition to ()2.20|) . Hence 

d+ 1 < diniN^) < d+2. (2.21) 

This extra zero eigenfunction is spherically symmetric and is also a zero eigen- 
function of the ordinary differential operator 

A M)0 = -A r + /i-/'(^) (2.22) 

on L 2 (IR + , r d ~ l dr), where A r is the radial Laplacian, 

A r = 8 2 r + d —^d r . (2.23) 

r 

Thus 

N(£,J = span{(0, ^), (<9 2 .^, 0), (N(^, ), 0)}. (2.24) 
For local nonlinearities ()2.10|) . N(A At ) — {0}, if either d = 1 or 

h'(r) + h"{r)r 2 > 0. (2.25) 

Alternative conditions on h for d > 1 are given in Appendix O Thus, for local 
nonlinearities, if either d — 1 or ()2.25j) holds, then condition |(F)| is satisfied (see 

also ISOl HZ1 E2I)- 

In any case this extra degeneracy, if it happens, is non-generic. If ()2.2|) 
holds for some /, then it also holds for small perturbations of /. On the other 
hand we expect that if (j2.2j) fails for a given / then there are arbitrarily small 
perturbations of / for which (|2.2|) is satisfied. 
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It is easy to check that Condition |(A)f|(Fj| are satisfied for nonlinearity (|1.2j) . 

Let S:=l d xl d x [0, 2tt) x / and S := R d xR J x [0, 2n) x J , where I is 
a closed, bounded interval on the positive real axis, specified in Condition |(C)[ 
and Jo is a closed interval contained in I\dl. Recall e := e v + e . Our main 
result is 

Theorem 2. Assume (jl.llj) awd |(04^[|(7^| are satisfied. Given e > and e v > 

defined in smc/j #iai e C 1 and an initial condition i[) satisfying 

Ue-W^o-T&Jlta < e o, (2.26) 

/or some o~o := {a , t>o, 7o, /-to} £ ^o- T7ien tdere zs a constant T > (indepen- 
dent of e v ,e but possibly dependent on I Q ) such that for times t(e v + e 2 ) < T 
the solution to with this initial condition can be written in the form 

ij;{x, t) = e l (^-(*-«)+7) - a )+w(x-a,t)) (2.27) 

w/iere 

Hk = 0(e) (2.28) 
and t/ie parameters v, a, 7 and /i satisfy the differential equations 

^ = -(W)(a) + 0(e 2 ), (2-29) 
d = t» + C(e 2 ), (2.30) 

^ = yU -^(a) + ^ 2 + 0(e 2 ), (2.31) 
/t = C(e 2 ). (2.32) 

The proof of this theorem is given in Section El 



3 The Hamiltonian and variational structure of 
generalized nonlinear Schrodinger equations 

It is well known that the generalized nonlinear Schrodinger equation is a Hamil- 
tonian equation of motion on an infinite-dimensional phase space. In this paper 
we make extensive use of this fact and of the symplectic geometry of phase space, 
in particular of certain "submanifolds" in this space. We outline briefly some 
facts which are relevant for us. 
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In this section we consider the generalized nonlinear Schrodinger equation 
eq. (HD), 

idfiP = (-A + V)iP - , 

under Conditions |(A)| and |(B)| only. 

We study eq. on the space Hi(IR d ,C). This space, considered as a real 
space (Hi(R d , R 2 ) = Hi(R d , R) ©Hi(R d ,R), ip ^ (Re ip, Imip)), and equipped 
with the symplectic form 

u>(u, v) := Im J uvd d x (3.1) 

(defined for u, v in the tangent space at a given point in Hi(R d ; M 2 ), which we 
can identify with H^IR^M 2 ) in our case), is a symplectic space. 
Define the Hamiltonian functional on Hi(M d , C) as 

HvW) -=\f | W| a + V\1>\ 2 d d x - Fty), (3.2) 

where F(i/j) is as in Condition |(A)[ i.e., F'(ip) = f(^p)- With these definitions 
eq. (jl.lj) can be written as 

= m'vW, (3.3) 

where J is the operator on H 1 (M d ,M 2 ) (strictly speaking J : (T^Hx)* i— > T^Hx) 
given by 

in block-diagonal notation. Thus J is a complex structure on Hi(IR d , R 2 ) corre- 
sponding to the operator i _1 on Hi(M d ,C). 

Observe that the space H^IB^, C) also has a real inner product (Riemannian 
metric) 

(u,v) := Re / uvd d x, (3.5) 

so that oj{u, v) = (u, J~ 1 v). 

Since the Hamiltonian H,y(ijj) is autonomous (it does not explicitly depend 
on time t) and invariant under the gauge transformation Hy(e' 7 ^) = 7iy(ip) for 
all 7 G [0,27r), we have conservation of energy, 7iy (■?/>) = const, and "mass" or 
"number of particles", = const, where 

J*W)'=\ J (3.6) 

12 



under the evolution of eq. (jl.ip . 

Due to Condition |(D)| the solitary wave profiles rj^ described in Condition 
are local minimizers of the Hamiltonian 7iv=o{ip) restricted to the spheres 



{ij G Hi : N(i>) = m} (3.7) 
for m > (see [20], Thm 3). Hence they are critical points of the functional 

£M'=\ j \^\ 2 + ^\ 2 d d x-F^), (3.8) 

where \i = n{m) is a Lagrange multiplier, and (jl.5j) is just the Euler-Lagrange 
equation for S^(ip). 

Observe that the functional also arises as (ti — t o )^(0) = 5y =o (0e 1M *) 

for any 4>(x), where Sy is the action for eq. fjl.ljl : 

SvM ■= I' (\ [ Imip$d d x + H v (ip))dt. (3.9) 

J t jRd J 

The functional S^ip) will play an important role in our approach. We will 
use it as a Lyapunov functional in estimating the fluctuations w. Finally, we 
note that the operator C v that appears in Condition |(F)| is the Hessian of 
at Tjn'. C v := S'^). 

4 Symmetries, zero modes, and the manifold of 
solitary waves 

In this section we introduce the manifold of solitary waves which is obtained 
by applying the generalized symmetry transforms (see below) to a fixed solitary 
wave. The tangent space to the manifold is introduced, and its inherited sym- 
plectic form is derived. Furthermore we prove the key fact that the inherited 
symplectic form is non-degenerate. 

Starting from this section we will often use the abbreviation 77 = 77^. 

Eq. (jl.lj) with V = is invariant under spatial translations 7^ tr , gauge trans- 
formations Xf, and Galilean transformations 7^ s , where 

T Q tr : i[>(x, t) ^ i>(x -a,t), T* : i/;(x, t) ^ e^(x, t) , (4.1) 

T/ al : i;(x, t) i-> e^-^^ix - vt, t) . (4.2) 
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(Transformations (|4.1j) - (|4.2|) map solutions of eq. with V = into solutions 
of (jl.ljl with V = 0.) The symmetries 7^ tr , Xf, T& al , have associated conserved 
quantities: field momentum, mass and 'center of mass motion', (cf. |17j). 

uPVtpd d x, J\ip\ 2 d d x, J Tp{x+2itV)ipd d x. (4.3) 

For f(ip) = \ip\ 2s tp and V = 0, eq. (jl.lj) is also invariant under the scaling 
transformation 

T; : tfj(x, t) ^ ^(y/jix, /it) . (4.4) 

When the external potential is introduced into the problem, the translational 
and Galilean invariance are broken. In particular, conservation of the field mo- 
mentum is replaced by the following 'Newton's law' (Ehrenfest's theorem) 

«W,-iV^) = -^,(W)^>, (4.5) 

which plays an important role in our analysis, and which is proved in Ap- 
pendix |XJ 

Let 7^ b : ip(x) i— > e2 v ' x ip(x), be the boost transform. We introduce the 
combined symmetry transformations S av ^: 

p(x) ^ p avi := S avy p = T^Zfpix) = e^ v < x - a ^p(x - a) . (4.6) 

Let T] avfl7 := S avi rj^. We define the manifold of solitary waves as 

M s := {rjavw ■ a, v, 7 , /i G R d x R d x [0, 2vr] x /} . (4.7) 

The tangent space to this manifold at the solitary wave profile 7/ M G M s is given 
by 

T^M S = span(z t , z g , z h , z s ) , (4.8) 

where 



z t :=V a r a tr ^| a=0 = -V^, (4.9) 

, (4.10) 
z h :=2V v T v ^ l r ]fl \ v=0t=0 = ix % , (4.11) 



9 <rs 



7=0 

al„ 
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and 

z s := . (4.12) 

Symmetries of (jl.lj) which are broken by the solitary wave solution 77 M e 1AI * 
lead to zero modes and "symplectically associated zero modes" of the Hessian 
C v . Namely we have for rj = 77^ 

C^Zt = , and C v z g = , (4-13) 

and 

£ v Zb = 2iz t , and C v z s = iz g . (4-14) 
The functions z h and z s are zero modes for the operator (i£,,) 2 , i.e., 

i£ v (i£ v Zh) = , and i£ v (i£ v z s ) = . (4-15) 

The relations (J4.13|) are proved by taking the derivatives of the equation 
E'^iT^T^rj^) = with respect to the parameters a and 7, at a = and 7 = 0, 
and similarly for ()4.14|) . 

We have shown above that 

T„M S C N gen (J/g , (4.16) 

where 

N gcn (JA,) := span{|J N((J£,)")} . (4.17) 

n>l 

In what follows we denote o := {a, f,7,/x} and r\ a := r] avifl . The manifold 
M s inherits a symplectic form from (Hi,uj). This symplectic form is determined 
by the operator 

n£ := h v „M = V n J~ X V na \t v „m s , (4.18) 

where V rja : T^Hi 1— > T^Mg is the L 2 -orthogonal projection onto T^ CT M S . 
Namely u> v<T (u, v) := (u, Qz^v). The key fact here is that this symplectic form is 
non-degenerate, i.e., the operator J -1 It^Ms, is invertible Wr] a e M s , as shown 
in Lemma ^ below. Define 



m(n) := 




(4.19) 



Lemma 1. Ifm'(fi) > 0, then VL^ is invertible. 
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Note that m\n) > is exactly the assumption for stability in Section |2l 



condition (TJ 



(4.20) 



Proof. We prove that Q v 1 is invertible by showing that its matrix D, 1 1 

((zj, J^ 1 Zk)) in the basis {zi, . . . , ^+2} = {zt, Z\>, z g, z s} in T^Mg is invertible. 
We compute 

/ -mil \ 
mil 
w! 
\ -w! J 

This matrix is invertible provided m' 7^ (m > always). Since f2~ CT x is related 
to Q^ 1 by a similarity transform (see (|5.6|l ). it is invertible as well. □ □ 

Corollary 2.1. £ T^Mg, there exists at least one element, z' £ T^Mg, sitcA 
i/iat uj v (z, z') 7^ 0. 

Proof. This follows directly from the non- degeneracy of £t~ l - □ □ 

Remark 1. The real form of the vectors ()4.9|) - (j4.12J) is 

(-V77,0), (0,7/), (O.aaj), («V/,0). (4.21) 

VKe abuse notation and denote both the real and complex representations of the 
vectors as —Vr], — Jr\, —Jxr/ and d^r], where for the real representation we 
interpret rj as (77, 0), and J is of the form ()3.4|) . whereas in the complex notation 

J = r\ 

Remark 2. For the special case with f{ip) = X\i/j\ 2s i/j, 

z s := d,T; V \^ =l t=Q = ~(~ + x ■ V)r? = . (4.22) 

Remark 3. Equation (ll.ljl with V = and wi/i a nonlinearity of the form 
()2.12|) . where W is spherically symmetric has also rotational symmetry: ip(x, t) \— > 
ip(Rx,t), R £ SO(d). T/izs symmetry does not play a role in our analysis since 
we consider only spherically symmetric solitary wave solutions of (jl.5|) . 

Remark 4. Note that for f(ip) = A|^| 2s, ?/ ; we can put the \x dependence on the 
same footing as the transformation induced for v, a and'-f. Indeed, define 

Sair/n = Sawy ° (4.23) 

Then rj av ^^ S av ^^rj\ . 
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5 Skew orthogonal decomposition 

In this section we introduce the skew orthogonal decomposition of a solution ip 
along the manifold M s and in the skew orthogonal direction (see [2]), and show 
that this uniquely defines the solitary wave solution parameters appearing in the 
decomposition. 

We recall the notation a = {7, a, v,/i} so that r\ a := r] aviiJi . Define the 5- 
neighborhood 

Us = W> G Hx : inf \\iP - r} a \\ Hl < 6} (5.1) 

of the manifold Mg := {rj a : a G So}, where 

S := M d x R d x [0, 2tt) x J . (5.2) 

Here Jo is any bounded closed interval contained in I\dl, (for the definition of 
/ Condition [fcj]). 

The main result of this section is 

Proposition 1. Let ip G Us- For 5 > sufficiently small, there exists a unique 
a = cr(^) G C 1 (?7 ( 5, S) such that 

uj(ij) - T] a ,z) = i.e., (ip - rj a , J~ x z) = , \/z G T^ CT M S . (5.3) 

Proof. We will use the following notation: 

cr := {j u a 2d +2} ■= {a, v, 7, //}, (5.4) 

,1; •••! z ^,2d+2 

h (5-5) 

and 

d k = d ak , k = d+ 1, ...,2d+ 2, ^ = ^ + ^^4.^^, fc=l, (5.6) 

Denote := S avi z^j so that (^.77^ = 2^ = S avi z^k- Clearly {z a j} is a basis 
in T,„\L.' 

We use the implicit function theorem for the map G : Hi x £ 1— > M M+2 , 
defined by 

G^, a) := (V - 7? CT , J-^ CTJ ) Vj = 1, 2d + 2. (5.7) 

We verify that (z) G G C 1 , (ii) G(r) ao ,ao) = for any er G £ and (m) 
5 CT G(r/ CT0 ,cr)| CT=CTo is invertible. 
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G is in C 1 in a, since so are rj a and z a j, and G is C 1 in ip since it is linear in 
ip. So (i) follows. 

(ii) follows directly from the definition of G. 
To prove (in), we use f!5.6|) to compute, 

d k Gj(Va , V)\ a=a0 = ~{Za ,k, J'^aod) = -<*«>,*> J^Z^k), (5.8) 

where in the last equality we used that S avi is symplectic (see Remark below). 
With the definition of f]" 1 in Section 0] we find 

d a G(r] ao ,a)\ = -ft" 1 ] 2d+2 . (5.9) 

Here we use the notation fi" 1 1 . , 2d+2 to denote the matrix of ft 1 in the basis 

{^ ,fc}fcS 2 - Thus d a G(rj mv a)\ a=(jQ is invertible for all a by Lemma [T] This 
shows (Hi). 

With properties (i)-(iii) the implicit function theorem implies that there 
exists a unique C 1 map o = cr(ip), satisfying G(if), cr(if))) = in a neighborhood, 
V CT0 , of r} ao . 

Now take a = {0,0,0,/i }, with fi Q G I. Denote n := {a, v,j} so that 
«S n := S av7 . Then for all n the map a, defined on iS^V CT0 , 

a(if>) :=noa(5 n -V) (5.10) 

where h o = {fi o n, fi} and n o n is defined by 5^ o S n =: iSnon, solves 

equation (|5.3jl for any -0 G 5nV CT0 . Since the neighborhood SnV ao with n G 
[0, 27r) and /io G / covers the neighborhood U$ the statement of the 
proposition follows. □ □ 

So if we know that for a given initial condition (jl.ljl has a C(R, Hi(M d )) fl 
C 1 (M,H_i(IR d )) solution ip = ip(t) which, for times < t < T, stays in the 
neighborhood Us, then by Proposition ^ the solitary wave solution parameters 
a trace out a unique C 1 trajectory a(ip(t)). We make the choice e <C 5. 

Remark 5. The operator S avi is a symplectic, or canonical, operator; i.e., it 
leaves (■, J -1 -), unchanged: (S avj w, J~ l S avi q) = (w,J~ 1 q). This follows from 
the fact that [J, £wy] = 0. Since z G T 5a ^ M M s implies that S~*z G T^Mg, we 
can write ()5.3|) in the form 

(S~^ - Vfl , J~ l z) = , Vz G T, M M S . (5.11) 

where r/^ and z depend on fi(t). 
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6 Equation of motion in the moving frame 

Given a solution ip to (jl.ljl . we define the parameterization {a,w} for it by the 
equations 

$ - v^J-^Ms and w := S' 1 ^ - Va ). (6.1) 

By Proposition Q this parameterization is well defined as long as ip G U$. In 
this section we find the equation for the parameters {a,w}. To this end we use 
the equation for the function u, defined by 

u := S- av \^. (6.2) 

We introduce the anti-self-adjoint generators 

Kj = d Xj , K d+j = - Jxj, JC 2d +i = -J, IC 2d +2 = d^, j = 1, ...,d (6.3) 

and the corresponding coefficients 

aj = hj - vj, a d+j = --Vj - d x .V(a), j = 1, d, (6.4) 



«2d+i = fi-^v 2 + ^a-v - V(a) - 7, a 2d+2 = -fi- (6.5) 



Let 



2d+2 2d+l 

a ■ JC = a jfcj an d a • JC = a j^j- (6.6) 
j=i i=i 
The main result in this section is 

Lemma 2. If if) satisfies (11.1)1 then u satisfies 

u = J((-A + //)u - f(u)) +a-]Cu + JK v u, (6.7) 

where 

TZ V := V(z + a) - V(a) - W(a) ■ x = 0((e v x) 2 ), (6.8) 

Proof. For this proof, we will use complex notation, i.e., J = i , and in partic- 
ular we note that [J, S avj ] = 0. Let if) a (x) = ip(x + a), and let <f) — \v ■ x + 7. 
With this notation u = e~ 1( ^/v Differentiating ()6.2|) . we find 

m = e~^(-A + /i + K)^a - f{ipa) + ie-^d • W>a + (^ ■ x + 7 - //)« , (6.9) 
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where V a (x) = V(x + a) and e 1( ^f(ipa) = f(u), by Condition |(B)| We rewrite 
the potential, using (|6.8j) . in the form 

V(x + a) = V(a) + VV(o) • x + TZ v {x). (6.10) 
= V(e-^a) + le-^aV^ (6.11) 

A(e-^a) + iV0 • V(e" 1 Va) + ie"^V0 • V^ a (6.12) 
A(e" 1 Va) + 2iV0 • V(e" 1 Va) - IV^V 1 *^ , (6.13) 

where = |t> -x + 7, to conclude ()6.7|) . Note that the parameters and their time 
derivatives, as well as V(a) and VV^(a) are collected into a. □ □ 

We now re-parameterize the non-linear Schrodinger equation into separate 
equations for a and u\ Recalling that u = r\ + w, d t r] = fid^r], £'^{rj) = 0, and 
S'^rj + w) = C v w + N v (w), where C v = S'^{rf), we can rewrite (|6.7|) as 

fid^T] + w = JC v w + JN v (w) + a ■ JC(r] + w) + JTlv{r] + w) (6.14) 

We collect the linear terms acting on w as C^ a w, where 

C v ,a := C v + TZ V + J- l a ■ IC, (6.15) 

and the remaining terms into a source term 

q(a) := a ■ JCi] + JlZyT], (6.16) 

to obtain 

w = JC v ,„w + N v (w) +q{a). (6.17) 

The equations for the parameters are obtained using the skew orthogonality 
condition and ()6.14|) . Let z G T^M,,. Upon recalling that (Jz, JC v w) = 0, and 
so 

= d t (Jz, w) = (Jz, w) + /i( Jd^z, w), (6.18) 

we find 

fi(Jz, d^r]) — a- (Jz, K/rf) = ^(Jd^z, w) + (z, N v (w))+a- (Jz, JCw) + (z, IZvijj+w)) . 



Now use 
and 
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Recall that d^r] = /C 2 d+2 ? 7, so the left-hand side is ^2j{Jz, K.jT))a,j = J2j(J z > z i) a h 
where we used KjT] = Zj. Now let z be one of the basis vectors, z = z k - Then 
the inner product on the left-hand side coincides with the definition of (fi" 1 )^-. 
Furthermore, using JC* = — /C, and [/C, J] = 0, we combine a ■ (Jz,JCw) and and 
^(Jd^z, w), into a ■ (JCz, Jw). The result is 

2d+2 

J^i^kjOj = (z k , N v (w) + n v {w + 7])) + a- (Kz k , Jw). (6.19) 

3=1 

Replacing a by the explicit expression, ()6.19|) reads, for k = 2d+l and k = 1, .., d, 
fi = {m'^y 1 ((t], JN v {w) + JK v w) - a ■ (JCy, w)) , (6.20) 
^v k = -d Xk V(a) + (m(/i)) _1 {{d k r), N v (w) + TZ v w) - a ■ (JCd k r], Jw) 

+(d kV ,K v ri)) , (6.21) 

where m(/i) := 2 _1 ||r]|| 2 , and we have used (Jr],7l v ri) = 0. For k = d + 1, ...,2d, 
k = 2d + 2, we use the expressions for fx and v k obtained above to find 

d k = v k + (m^))" 1 [(x k T], JN v (w) + JTl v w) + a ■ (ICx k r], w)} , (6.22) 
7 = - j^ 2 + \a • v - V(a) - (m'^))^ 1 [(d^r], N v (w) + TZ v w) 
-a ■ (JCd^T], Jw) + (d^jTlvr))] , 

where we used (xkr], JTZyf]) = 0, and rj(x) = rj(\x\) so that (xt],7]) = 0. Fur- 
thermore, observe that all terms containing w and IZy are of higher order. We 
abbreviate this as follows 

a = X(a) -5X{a,w). (6.23) 

For the estimates used later, we note that dj — Xj(a) = aj. We formalize the 
above calculation in the following Proposition: 

Proposition 2. (1) The parameters a and the fluctuation w (defined in (|6.1|) ) 
satisfy the equations 

& = X(a)-6X(a,w), (6.24) 

and 

w = JC Vta w + JN v (w) + q(a). (6.25) 
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Here C Vt(T , q(a) and N v (w) are given by ()6.15|) . ()6.16j) and (|2.5|) respectively, 
and (with defined in LemmaU}) 

2d+2 

5Xj(a,w) = 22{tt v )jk[(zk, N v (w) + TZ v (w + rf)) + a ■ (JCz k , Jw)), (6.26) 

k=l 

Vj = 1, 2d + 2, with {zk}ltti '■= { z t, z g , z s }, and IZy given by ()6.8|1 . 
(2) The vector field SX satisfies the following estimate for ||w||hi — 1- 

6X = 0{\a\\\w\\ + e v 2 + \\w\\ 2 Ul ) , (6.27) 

where \a\ := maXj = i v .. j2( i+2 |ci!j|. 

Proof. (1) is the result of the calculation done in ()6.14|) - (|6.23|) . In particular 
HHU) follows from 

(2) Estimate ()6.27|) follows directly from ([6.26)1 . together with the facts that 
\\K v z k \\ = 0(e v 2 ) and H^HUh-! < c||tw||^ for ||w|| Hl < 1 (see □ □ 

The goal is to show that sup tg ( 0T ) \SX\ = 0(e v 2 + e 2 ) and sup tg ( 0T ) HHIhi — 
0{e v + e ), for some T = 0{l/{e v + e 2 )). 

7 Approximate Conservation of a Lyapunov func- 
tional 

In this section, we show that the Lyapunov functional £^{u) — £^(77^), is approx- 
imately conserved. Recall r\ = rj^ is the solitary wave profile (see Section |TJ) and 
u is the solution ^ of eq. (jl.ljl . transformed to the moving frame: w := S~J^ip 
(«S a „ 7 is defined in Section HJ). We use the skew-orthogonal decomposition of w 
(see Sectional): u := S^Lif) = rj + w, with 

0^(11?, = for all z G T^M S , (7.1) 

provided ip EU$. The main result of this section is the following 

Proposition 3. Let ip G Us solve eq. and let u, w and 77 be defined as 

above. Then 

- £ M (r/ M )) = OflalHI^ + e^lMk + CvHIhJ- (7-2) 
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To prove this proposition we use the following 
Lemma 3. Let u be defined as above. Then 

dtS^u) = ~/i|M| 2 - ((-?) + VK)m, Vu). (7.3) 

of Proposition^ We first recall that 77 is a critical point to S^u), thus 

d t £M = ^||r/|| 2 . (7.4) 

Using Lemma El we find 

d t (S, l (u)-S^ V )) = A-B. (7.5) 

where A := 2 1 /x ( 1 1 1 1 2 — ||r/|| 2 ) and B := ((2 _1 ?) + VK)m, V«). First we use the 
decomposition u = 77 + w, the condition = (iz g ,w) = (t),w) (from 1)7. ljl ). and 
the estimate \fi\ — |a 2 d+2 1 < | or | to obtain 

A = ^fi\\w\\ 2 = 0(H|M| 2 ). (7.6) 

For the term B, recall that = (iz t ,w) = (iVr],w), and furthermore that 
(1577, Vr/) = for any real- valued function q £ L°°. Then 

B = ((±v + W a )iw,Vw) + ((VKH W) + ((VV a )w, Vr/). (7.7) 
Now, we use that VV(a) • (iu>, Vr/) = = W(a) • (177, Vto), to obtain 

B = (-v + W(a)) ■ (iw, Vw) + ((W a - W(a))iw, W)+ 

((W a - W(a))ir/, Vw) + ((W a - W(a))iw, Vr/). (7.8) 

Recall that — aj+j := + 9 3 -V(a) (see (16.5)1 ) . Since |a 3 -| < |a|, the first term on 
the right hand side is 0(\a\ IIHIhJ- Since W a = 0(e v ) the second term on the 
right-hand side is 0(ey ||?x>||j![ L. Finally due to W a - W(a) = C(e v 2 |x|) and 
|x|r/, |x|Vr/ G L 2 , (Condition |(C)| ) the third and fourth terms are ©(ev^HioHHj- 
Collecting these estimates we arrive at 

B = 0{\a\\\w\\ 2 Ui + e Y 2 \\w\\ Hl + ^NIhJ- (7-9) 
Relations (J73D, (J7SJ) and imply (|72J . □ □ 
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It remains to prove Lemma 01 To this end we note the following 
Lemma 4. Let ip be a solution to eq. (jl.lj) . Then 

dt\ f Vh/f=((VV>A,V^>. (7.10) 

Proof. ()7.10|) is obtained by integrating the relation 

d t (V\^\ 2 ) = iV ■ (V$Vip - VtpV^) - i(W) ■ (</>W - VW), (7.11) 
which follows from the nonlinear Schrodinger equation (jl.lj) . □ □ 

of Lemma\^ Note that the identity 

n v (S-\^) + ^||5-^|| 2 -\f V\S~ av \^ = £,{u). (7.12) 
holds for all u = S^Lip. We observe the following relations 

II^VII 2 = ||V>|| 2 , / V\S-^\ 2 = J V-M 2 (7.13) 

and 

2H V {S-^) = 2H v (iP) + ^ 2 ||^f - v ■ W> + / - VM\ 2 (7.14) 

Inserting the above relations into (j7.12jl we find 

2H v (i>) + (\v 2 + t i)W\\ 2 -v{ii>,V^)- j V\^\ 2 = 28,{u). (7.15) 

We now take the time derivative of the above relation. Using the fact that 
Tiv(i/)) and HV'II 2 are conserved quantities, and using LemmaEJand Ehrenfest's 
theorem (j4.5j) we find 

(h-v+fi)\\^\\ 2 -v-(iij, V^}+v(VV^,^}-2((VV)uP, V^> = 2d t S^u). (7.16) 

Collecting terms of the form (|?) + VV) and replacing ip with S avi u gives 

fi\\u\\ 2 + (i(t> + 2VK) • Vu, u) = 2d t E li {u). (7.17) 

By observing that (qiVu,u) = — (qvu, Vif) for real g G L°°, we arrive at the 
result of the lemma. □ □ 
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8 Lower bound on the Lyapunov functional 

Proposition 4. Let 77 and w be as in Proposition^ Then there exist constants 
p > and c > independent of e v and e such that for ||iu||hi < 1> 

1^(77 +w)- ^(77)| > ^\\w\\ 2 ni - c\\w\\^ . (8.1) 

Proof. We expand S^u) around 77. Using that 77 is a critical point to £ M (S'Jr)) = 
0), we write 

S,( V + w)- £M = ±(w, C v w) + R®(w) , (8.2) 

where, recall, C v := S'^ijf) and where R v (w) is defined in fl2.7|) . From condi- 
tion |( A) | and ()2~1) we have for ||tw|| Hl < 1 

\R®(w)\<c\\w\\ 3 Bl . (8.3) 

Let := {w G Ri(R d ) : ^(w,z) = 0,Vz G T^M S , ||u;|| Hl = 1}. It is shown 
in Appendix |D] (see [IB]) that 

p := inf (w,C v w) > 0. (8.4) 

Hence, for w that satisfy 1)7.1 Jl we have the following coercivity estimate 

(w,£ v w) > pWwW^. (8.5) 

Using this estimate and the bound ()8.3|) on R v in ()8.2|) we arrive at (|8. 1|) . □ □ 

Remark 6. Since cr eS s(£r)) = [/•*, 00), we /iaue £/ia£ p < /i. H^e expect that for a 
wide class of nonlinearities p > c/i /or some constant c > 0. 



9 Upper bound on ||w||hi an d proof of the main 
result 

In this section we prove the main result of the paper, by providing an upper 
bound on HwIIhv To achieve this, we use both the approximate conservation 
of Sfj_ and the lower bound on the Lyapunov functional. For vector functions 
s 1— > w(s) G Hi, and s t— > a(s) G M 2d+2 we introduce the norms |||w|||hi := 
su Pse[o,t] j| j| Hi and \a\oo := sup s<t |a(s)|. We state the main result of this 
section 
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Proposition 5. Assume e v , e are sufficiently small. There are constants c, d < 
oo, independent of e v and e such that for t < c(e v + e 2 ) -1 , 

II^IIhi < c'(e v + e ), (9.1) 
|«|oo <c'(e v 2 + e 2 ), (9.2) 

where w = S^Lip — rj^ and aij = &j — Xj(o~) (see (|6.24j) ). 
Denote 

A£ := £^{u Q ) - £ M0 (?7 M0 ) , (9.3) 
where u := S^ift | . We begin with two simple auxiliary lemmas. 

Lemma 5. There exists a constant c > independent of e v and e such that 
\A£\ satisfies the inequality 

\A£\ < ce 2 . (9.4) 
Proof. To estimate A£, we use the fact that £„ (^ Mo ) = to obtain 

2A£ = 2(£ IM) (r) tl0 + w ) - £ Mo (^ ) - (Z'^iVw), Wq)) 

= ||Vwo|| 2 + /loINf - #f (9-5) 

Since ||wo||hi < 1, the estimate ()2.4|) gives 

|l4 2) WII<c|ho|lL. (9.6) 

This together with the identity (|9.5|) and the fact that ||u>o||hi < e o < 1 gives 
the estimate JOj)- □ □ 

Lemma 6. Let p > be the coercivity constant given in ()8.4|) . There exists a 
constant c independent of e v and e , such that for ||w||hi < 1 

HIMIIhi < ce o 2 + ct(e^ 2 1 1| || | Hl + (ev + lalocOHMIlHi) + c HMIIhi ( 9 - 7 ) 

Proof. Using the approximate conservation of the Lyapunov functional (the time 
integral of ()7.2j0 . and pulling sup of norms out of the time integral we obtain 

\£^ + w)-£M\ < \A£\ + ct(e v 2 \\\w\\\ Hl + (e v + laUIHHiiJ ( 9 - 8 ) 

where \A£\ is defined in (19 .3|) . Now we substitute ()9.8|) into the lower bound 
(|8.1|) . use the initial condition estimate ()9.4|) and take the sup se j 0jt ] of the result- 
ing expression to obtain (|9.7j) . □ □ 
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Figure 1: Schematic graph of the left-hand side of (|9.9[) . 

of Proposition Using the triangle inequality we derive from (|9.7|) that the 
function X{t) := IH^HIfi! satisfies the equation 

< pe v 2 + ce 2 - pX 2 + cX 3 (9.9) 

for times ct(e v + |a|) < p/4, provided X < 1. The graph of the right-hand 
side of f!9.9|) is shown in Figure [T] Thus for pe v 2 + ce 2 < q , see Figure [U 
X = IHHIIhi < P if ^lt=o — Pi where p is the smallest positive zero of the 
left-hand side of ()9.9|) (see Figure HJ), provided p < 1. For e + e v sufficiently 
small, c'(e v + p _1 / 2 e ) < p < c(e v + p _1 / 2 e ) and therefore 

IMIhi < HMIIhi < c{e v + — e ) (9.10) 

provided HwoHhj < 1- Substituting (j9.10j) into ()6.27|) we obtain from ()6.24|) and 
the notation a := a — X(a) that 

laU <c( ev . 2 + p-V), (9.11) 

for times t < c(e v + p _1 e 2 ) -1 . This completes the proof. □ □ 

A Ehrenfest's theorem 

In this appendix we prove Ehrenfest's theorem, eq. ()4.5|) . Denote I^{t) := 
(ifi,-id Xj il)). Let V solve eq. (HHJ and assume that i[) G C(R,H 2 ) H C^M.L 2 ). 
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Then I^(t) is in C 1 , and after a straightforward calculation using eq. we 
find 

4 = / -|^|\yd d x =: L^t). (A.l) 



By the fundamental theorem of calculus, we obtain 

I4t) = l4t )+ [ Lip{s)ds. (A.2) 

Now take an initial condition ip G Hi. Pick ^ ,n £ H 2 such that ipo,n — ► 4>o 
in Hi. Then by Theorems 5.2 and 4.2 of [TT] . the solutions ^ n corresponding 
to the initial conditions ^ ,n satisfy ^ G C(R, H 2 ) PI C 1 (R, L 2 ) and ip n — > ip in 
C(R,Hi). Thus since 

hAt) = hM+ [ L^ds, (A.3) 



we have 



I f (t) = I f (t )+ [ L^ds. (A A) 

We furthermore observe that both I^(t) and L^(t) are continuous in t for ^ G 
C(R,Hi)nC 1 (R,H_i). Hence flA~!) implies that i^(t) is C 1 and satisfies /^(t) = 
L f (t). This proves □ 



B Minimization under constraint and spectrum 
of Hessian 

In this appendix we show that the operator C n := S'^(r]) has exactly one negative 
eigenvalue. The argument below is probably well-known, but we did not find 
it in the literature. Let X be a Banach space and K G C 3 (X, R) be a given 
functional. Define the set 

M = {u G X : K(u) = 0} . (B.l) 

We have the following 

Proposition 6. Let £ be a C 2 functional on X. Assume there is a Hilbert 
space, H, such that H D X ; densely, and that the Hessian quadratic form 
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Hess £(u)(a, (5), a, (3 G X ; defines a self-adjoint operator £"{u) on H such that 
(a,£"(u)(3) = ttess£(u)(a,/3), Va 6 X J 6 D(£"(u)) c x - Let V be a mini- 
mizer of £ on the set M C X defined in flB.l|) . Assume K'(j]) ^ 0. T/ien the 
Hessian operator £" '(rj) has at most one negative eigenvalue. 

Proof. Let rj be a minimizer of £ on M. Then rj satisfies 

£' (77) = and £"(77) > . (B.2) 

Here S'(ij) : T V M ^ R and £"{rj) : T V M ^ (T,M)*, where 

T„M = {9 s r7 s | s=0 : r/ s G C^fO.e^M), r/ s=0 = v} • (B.3) 

We claim that T^M can be written as 

T„M = {£ G X : (if '(77), = 0} =: tf'fo)- 1 . (B.4) 

Indeed, if £ £ T^M then there exists r/ s as in (jB.3|) with £ = d s \ s=Q r] s and 
therefore 

= ^(77 s )| s=0 =(if'(r7),O. (B.5) 

On the other hand if £ G K'irj) 1 - then we can find 7/ s such that r] s= o = r), 
9 s \ s=0 Vs — £ an d K{j] s ) = by solving the equation f(a, s) — where 

/(a, s) := \K( V + s£ + s 2 aK'{r,)) (B.6) 

for a and setting ?7 S = r)+s£+s 2 aK'(r]). The latter equation has a unique solution 
for s sufficiently small since f(b,0) = 0, where b = — (£,K"(t))£)/\\K'(ji)\\ 2 and 
d a f(b,0) = \\K'(rj) || 2 , and therefore the implicit function theorem is applicable. 
Now, the second equation in (jB.2|) can be rewritten as 



inf (t,e»(rj)$>0. (B.7) 

£6K'(?7)-L 

According to the max-min principle, the number of non-positive eigenvalues of 
£"{r]) is less or equal to the co-dimension of K' (rj) , which is 1. □ □ 

In our case X = Hi(R), H = L 2 (R d ), £ := £„, where 

£„{u) = i I \Vu\ 2 + fi\u\ 2 d d x - F(V>) , (B.8) 
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and rj is a minimizer to with the constraint 

K(u) -=\J H 2 d d x - m = . (B.9) 

This implies that C v := £^{t]) has at most one negative eigenvalue. 

Proposition 7. Under condition \(D)\ C V has exactly one negative eigenvalue. 

Proof. By Proposition |U] has at most one negative eigenvalue. One the other 
hand, since C^d^rj = —77, we have 

(dtf, C^rj) = -\dp( v 2 d d x < 0, (B.10) 



by Condition |(D)| Therefore by a variational principle C v has at least one 
negative eigenvalue. Thus C v has exactly one negative eigenvalue. □ □ 

C The null space of C v 

In this appendix we discuss condition |(F)| (see (|2.2|) ). We represent complex 
functions u(x) = u\{x) + iU2(x) as real vectors (ui(x),U2(x)). In this represen- 
tation the operator C v takes the form 

A = ( L ' r ) > (CI) 



L 2 

(on L 2 (R d , R) © L 2 (R d , R)), where 

L 1 = -A + /i-/%), (C.2) 

and 

L 2 = -A + f i-fV>{rj), (C.3) 

with / (1) (?7) := <9 ReV ,(Re /)(?y) and / (2) (r?) := 9[ m ^(Im/)(^). The diagonal form 
of C v follows from the diagonal form of f'{rj): 

f'{v) = f®(v))- (C.4) 

The latter follows from the relation f(T c ip) = T c f(ip), where T c is complex 
conjugation, and the fact that rj is real. This relation, in turn, follows from 
F(T c i/j) = F(ip) (see Condition p)]). 
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The matrix operator (f(TT| is then extended to L 2 (R d , C) © L 2 (R d , C). The 
operators L\ and L 2 are self-adjoint, with essential spectra given by 



0*a(Lj) = \pt,<x), j = 1,2. 



(C.5) 



Relation ()2.20j) implies that 



0^.77 e N(L0, Vj = l,...,d, 



(C.6) 



and 



v e N(L 2 ). 



(C.7) 



From now on, we assume that / is a local nonlinearity. The fact that rj > 0, 
implies by a Perron- Frobenius argument (sec |36j) that 



Thus it remains to analyze the operator L\. To begin with, we observe that 
since C n has exactly one negative eigenvalue, relation (jC.8|) implies that L\ has 
exactly one negative eigenvalue (the same as C v ). 

First we consider the case d — 1. Then the zero mode, 7/, of Li has ex- 
actly one zero (at x — 0) and consequently (by Sturm-Liouville theory) has 
exactly one negative eigenvalue, as we concluded above from general considera- 
tions. (Remember that the lowest eigenvalue in our case is simple, and that the 
nonlinearity is local) 

Theorem 3 ( 07]). Let d = 1. Then N(L X ) = Crj' . 

Proof. The proof follows Weinstein jlTj. We know one solution, rj'. The fact 
that i]"(0) ^ allow us to choose the first linearly independent solution to 
be w\ = r]'/r]"(0). Then w%(0) = 0, w[(0) = 1. Consider a second linearly 
independent solution, W2, with u^O) = 1. Since the Wronskian 



W[w 1 ,w 2 }(x = 0) = wi(0K(0) -w[(Q)w 2 (0) = -102(0) = -1 , (CIO) 



L 2 > and N(L 2 ) = Cry. 



(C.8) 



(C.9) 



is constant with respect to x, and since 



we have W / [tfi,w 2 ] = — 1 for all x, and therefore 



T)'(x)w' 2 (x) — 7]"(x)w 2 (x) 



V'(0) > . 



(C.ll) 
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The last equation can be rewritten as 
Now, for e > and x > 0, 



w 2 (x) = w 2 (e)^ - rj'(x)rj"(0) jf dr. (C.13) 



Since t/(x) < for x > 0, r]"{0) < 0, and \rf(x)\ < Ce~ c ^, we have 

w 2 {x) < w 2 (e)^ - i/(s)i/'(0) jf ^ -oo as * ^ oo . (C.14) 

Hence w 2 ^ N(Li) and we are done. □ □ 

Next, we consider the case d > 2. We use the assumption that ?y is spheri- 
cally symmetric. In this case the operator L\ (and also L 2 and C v ) is spherically 
symmetric; i.e., it commutes with the action of the rotation group SO(d) on 
L 2 (IR d ). As a result it can be decomposed into a direct sum of ordinary differ- 
ential operators corresponding to the eigenfunction expansion of the Laplacian 
Agd-i on the sphere S d_1 . Denote the orthonormal eigenf unctions of — Agd-i 
corresponding to the eigenvalues X k = k(d — 2 + k) by Y k {6), 9 e S d_1 . Then the 
operator L\ can we written as the direct sum 

Li = ®A^ k (C.15) 

acting on the direct sum L 2 (R d , d d x) = 0^ o L 2 (M+, r d_1 dr) <g> Y k (9), where the 
operators A^ k are defined on L 2 (R + , r d_1 dr) by 

A^ k = -A r + fi + V k (r), (C.16) 

and where A r = <9 2 + (d — l)r~ 1 d r is the radial Laplacian in M. d and V k (r) := 
-f'(r))(r) + A fc r~ 2 . Clearly A*^ k = A^ k and a css (A^ k ) = oo). 
Now observe that 

(d x .ri)(x)=x j rf(r) E L 2 (K + , r^dr) ® ^(0), (C.17) 

where x = e S d_1 . Hence 

V G N^x). (C.18) 
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Since rf (r) < for r > 0, we have by an extension of the Perron- Frobenius (or 
Sturm-Liouville) theory (see [SI]) that zero is the lowest eigenvalue of A^i and 
is simple. Thus 

A^x > and N(A M ) = Ct/. (C.19) 

Next, since for k > 2 

A^k ~ = > 0, (C.20) 

we conclude that A^ > for k > 2. 

Theorem 4. N(y4 Atj0 ) zs trivial provided f(ip) = h(\ip\ 2 )i/j with h satisfying 

h'(r) + h"(r)r > 0, r > 0. (C.21) 

Proof. Assume there exists ( G L 2 such that A^^C, = 0. Then is the second 
eigenvalue of A^ Q , and so the corresponding eigenfunction £ has exactly one zero 
in (0, oo), say at r . Observe the following properties 

i. £±Ran(v4 Atj0 ); 

ii. A^qT] = 2h'(i] 2 )r] 3 ; 
iii. A M)0 <9 M ?7 = -77. 

Properties Ell and Irnl follow from the equation 

(-A r + f i)r ] -f(r ] ) = (C.22) 

for 77. 

Properties RHml imply the relation 

(6, (h'(v 2 W - "to) = 0, Va 6 R. (C.23) 

Since rj is monotonically decreasing from some 77(0) > at to as r — > 00 
(see e.g., [I]), and since h'(s)s is monotonically increasing function of s by con- 
dition (|C.21j) . we can choose a such that the monotonically decreasing function 
h'((r](r)) 2 )r] 2 (r) — a has a zero exactly at tq. In that case the left hand side of 
flC.23|) is non-zero, which leads to a contradiction. Thus the equation A^o£ = 
has no nontrivial solutions in L 2 . □ □ 
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Finally we present the conditions under which McLeod's [2.9.^ uniqueness 
proof of positive solitons implies that A^ has trivial null space in the case 
d > 1. There exists a > such that 

ps — sh(s 2 ) > 0, for < s < a, 

ps - sh(s 2 ) < 0, for a < s < oo, (C.24) 
(sh(s 2 ) — [is)' > 0, when s = a, 

and for each S > a, 3A = X(S) G C((a, oo), R + ) such that 

K(s,X) > for s e (0,S), K(s,X) < for s G (S,oo), (C.25) 

where 

K(s, X)=/is + Xs 3 h'(s 2 ) - sh(s 2 ). 

D Coercivity of C v 

The goal of this appendix is to prove the following result, essentially due to |4*5] . 

Proposition 8. There is p' > such that if w satisfies u(w, z) = Wz G T^Mg, 
then 

(w,jC v w) >p'\\wf Hl . (D.l) 

Proof. We break the proof into three steps. The proof utilizes the fact that C rj 
has exactly one non-degenerate negative eigenvalue and the assumption |(F)| i.e., 
that N(£^) = span{(0, 77), (d Xj r],0), j = l,...,d}. 

Step 1. Let Xi = {w G H a : ||w|| = 1 , ((r),0),w) = 0}. T/ien 

inf (w, £„w) = . (D.2) 

Proof. Let a := inf^gXi (w, C n w), (see ()C.ljl ). Clearly v < a < 0, where z/ < 
is the negative eigenvalue of C v (see Appendix [BJ). That a < is clear, as 
w = (0,T])/\\t]\\ G Xi yields (w,£ v w) = 0. Moreover a ^ v. Indeed if a = v 
then the minimizer, v, of (jD.2|) would be an eigenfunction of C v corresponding 
to the smallest eigenvalue, v. Since (?7,0)_Lf, (77, 0)_LN(£, ? ) and since is the 
only negative eigenvalue of (see Proposition EJ), we conclude that (77, 0) is in 
the spectral subspace of corresponding to the interval [5, 00) for some 5 > 0. 
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Therefore C r] 1 (r], 0) is well defined and ((77, 0), C l {r] : 0)) > 0. On the other hand 
the equation £^(8^,0) = — (77, 0) implies that 

((r ]j 0),£- 1 (r ] ,0)} = -m'( f i)<0 (D.3) 

by Condition |(D)| which contradicts ((r/, 0), £^(77, 0)) > 0. Hence a = v is 
impossible. 

To show that a = we use the Euler-Lagrange equations corresponding to 
(El 

C v w = aw + f3(r],0) (D.4) 

where a and (3 are the Lagrange multipliers corresponding to = 1 and 
((r),Q),w) = respectively. Assume v < a < 0. If (3 = 0, then a would 
be a negative eigenvalue in (1/, 0) which contradicts that v is the only negative 
eigenvalue. Thus (3 7^ 0. Given v < a < 0, (3 7^ 0, we can solve the Euler- 
Lagrange equation as 

w = (3{C V - a) -1 (77,0). (D.5) 

The inner product of the equation above with (77, 0), the orthogonality relation 
(w, (77, 0)) = 0, and (3 7^ 0, give 

= ((77, 0), (£„ + |a|)- 1 (7 7 , 0)) =: ?(|a|). (D.6) 

q(X) is analytic in A G (0, |z/|), and hence different iable. Moreover, it is mono- 
tonically decreasing, since 

q'(X) = -((77, 0), (£, + A) -2 (77, 0)) = -\\(C V + A) _1 (77, 0)|| 2 < 0. (D.7) 

Furthermore by flE| g(0) = ((77, 0), £^(77, 0)) < 0. Thus q(\a\) ^ 0, for 
a G (^,0), which contradicts ()D.6|) . Hence a = 0. □ □ 

Step 2. Lei X := {w G Hi(R d , C) : ||w|| = 1, u{w, z) = 0, \/z G T^MJ. T/ien 

inf (w, £„w) > . (D.8) 

Proof. The Euler-Lagrange equation corresponding to (jD.8|) is 

Z^w = aw + 7 fc Jz fc (D.9) 
k 

where {z^} is a basis for T^Mg. Here a and {7^} are the Lagrange multipliers 
corresponding to the constraints \\w\\ = 1 and u(w,Zk) = V7c respectively. 
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Note that a = (w,C v w), and that X C X 1; hence a > 0. Assume that a = 0, 
and that one 7^ 7^ 0. Then for at least one z k G T^M S , we have 

(z k , C v w) = 7i(^r,)jfc 7^ 0, (D.10) 

for some k, which contradicts (z k , C q w) = (C v z k , w) = Wk. Here we have 
used that det Q v ^ 0, and that z k is either a zero-eigenfunction or an associated 
zero-mode for C v . Thus either a > or jj = 0. Consider the latter case. In this 
case 

£ v w = 0. (D.ll) 

which implies that w G N(£ J? ). Since N(£^) C T^M S , the relation u{w,z k ) = 
for all z k G T^M S , contradicts the non-degeneracy of (see Corollary 12. lj) . 
Thus a > 0. □ □ 

Step 3. End of Proof. Eq. (|D.8|) implies that there exists a p" > such that 

{w,£ v w) > p"\\w\\ 2 , (D.12) 

for some p" = p"{p). To improve the coercivity from L 2 to Hi, we let < 5 < 1, 
and estimate (w, C v w) using (|D.12J) as 

(1 - <5)p"|M| 2 + S{w, C v w) < (w, C v w) . (D.13) 

Upon using the explicit form of C v we find that 

(w, C v w) > \\Vw\\ 2 - CJw\\ 2 , (D.14) 

where 

C, = su P (p+\f'( V )\). (D.15) 

X 

The last two estimates with 5 := p"(l + p" + C M ) _1 imply 

(w,£ v w) > p'||w||hi , (D.16) 
where p' = p"(l + p" + C A1 ) _1 . This concludes the proof of Proposition |H1 □ □ 
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